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Abstract: 

It is proved that the Laurent expansion of the following Gauss hypergeometric functions, 



2 Fi (ii + ae, I 2 + be; I 3 + ce; z) , 
2 Fi (h + ae, I 2 + be;h + \ + ce; z) , 



><! ■ 

y ■ 2-^1 (h + i + h + be; I 3 + ce; z) , 

iF\ (Ji + \ + ae, I 2 + be; I 3 + \ + ce; z) , 

2 Fi (Ji + \ + ae, I 2 + \ + be;h + \ + ce; z) , 

where Ii,I 2 , 13 are an arbitrary integer nonnegative numbers, a, b, c are an arbitrary num- 
bers and e is an arbitrary small parameters, are expressible in terms of the harmonic 
poly logarithms of Remiddi and Vermaseren with polynomial coefficients. An efficient algo- 
rithm for the calculation of the higher-order coefficients of Laurent expansion is constructed. 
Some particular cases of Gauss hypergeometric functions are also discussed. 



Keywords: 


Gauss hypergeometric functions, harmonic polylogarithms, colour 


poly logarithms, Laurent expansion of Gauss hypergeometric function, multiloop] 


calculations 





'Supported by NATO Grant PST.CLG.980342 and DOE grant DE-FG02-05ER41399 



Contents 



[T]. Introduction |l| 

|2|. All-order e-expansion ||| 

|2.1| Non-zero values of the e-dependent part H 

|2.1.1| Integer values of of e-independent parameters [| 

|2.1.2| Half-integer values of of e-independent parameters || 



|A| . The Laurent expansion of Gauss hypergeometric functions with half- 



2.2 Zero-values of the e-dependent part of upper parameters 10 



|3[ Some particular cases 11 



A\ The generalized log-sine functions and their generalization 11 



7A Special cases: all-order e-expansion in terms of Nielsen polylogarithms 13 



Conclusions 14 



integer values of parameters around z = 1 15 



1. Introduction 

One of the most powerful techniques for calculating Feynman diagrams is based on their 
presentation in terms of hypergeometric functions. We will call this the hypergeometric 
function representation of Feynman diagrams. Such a representation can be used for nu- 
merical evaluation, construction of the asymptotic expansion, etc. One of the unsolved 
problems in this program is obtaining the proper representation for a diagram with an ar- 
bitrary number of legs and loops. Direct use of a- or Feynman parameters representations 
|l[ is not very helpful in solving this problem. The Mellin-Barnes technique is restricted to 
several topologies |2|, ||| [|, [j| . The negative dimension approach || has a similar restriction 
@. The most investigated diagrams are the master integrals (typically, integrals with the 
power of each propagator equal to unity). The differential Q and/or difference equation 
B techniques are usually used to obtain such representations. The known cases include 



the one-loop diagrams |1C| 11 1, two-loop propagator- type diagrams with special mass and 
momentum values fl2| , several three-loop bubble-type diagrams ||, three-loop vertex-type 
diagrams [13], and four- loop bubble- type diagrams |l4j]. For practical application how- 



ever, it is necessary to construct the of e-expansion (Laurent expansion) of hypergeometric 



functions. There is some evidence that the multiple polylogarithms |15|, [16], 17], 

£n\ z m 2 _ _ _ z m n 

Li/c 1 ,fe 2) ...,fe n (^1^2,--- ,Z n ) = Yl \ ^2 \- ( L1 ) 
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are sufficient for parametrizing the coefficients of the e-expansion of some, but not all 1 , 



hypergeometric functions [18]. 

In some particular cases, the result of the Laurent expansion can be written in terms 



of simpler functions. In particular, at the present moment, it is commonly accepted [19 



that the generalized hypergeometric functions with an arbitrary set of integer parameters 
can be presented in terms of harmonic polylogarithms pO]. The idea of the proof is based 



on the properties of nested sums [21]: the analytical coefficients of the e-expansion of any 
generalized hypergeometric function with integer parameters can be reduced to a set of 
some basic harmonic series of the type 



J2^S ai (j-l)---S ap (j-l), 
i=i 3 



r X2) 



where z is an arbitrary argument and S a (j) is an harmonic sum defined as S a (j) = Ylk=i W' 
Series of this type are expressible in terms of the Remiddi-Vermaseren harmonic polyloga- 
rithms. 

However for hypergeometric functions with half-integer values of parameters, the new 
type of sums, multiple (inverse) binomial sums [22, [|, ||, 23 1 are generated: 



J ai,- 



,a p ; bx, 



E 

3=1 



1 Z-? 
2j\ k f 



S ai U -I)-" S ap (j - l)S h (2j -!)••■ S bq (2j - 1) 



J1.3) 



For particular values of k, the sums (1.3) are called 



generalized harmonic 

1 inverse binomial 
-1 binomial 



sums . 



At the present moment, there is no proof that any multiple (inverse) binomial sums 
can be expressed in terms of harmonic polylogarithms only. This problem was investigated 
in Ref. p3] for multiple inverse binomial sums up to weight 4- In Ref. p4| , it was shown 
that some of the multiple inverse binomial sums are not expressible in terms of harmonic 
polylogarithms of simple argument. In Ref. [25|, the results of Ref. |23|] were extended 



on the case of special combinations of multiple binomial sums and multiple generalized 
harmonic sums. However, the Laurent expansion of a hypergeometric function in general 
contains combinations of multiple sums. These combinations may be expressed in terms 
of harmonic polylogarithms. From this point of view, the construction of the analytical 
coefficients of the e-expansion of hypergeometric functions can be done independently from 
existing results for each individual multiple sum. 2 



1 We are thankful to S. Weinzierl for this information. 

2 We are indebted to A. Davydychev for discussion on this subject. 
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The simplest hypergeometric function is the Gauss hypergeometric function 2 F\ (a, b; c; z) 
1 26, 27, p£j j. It satisfies the second-order differential equation 

Tz{4z +c - 1 ) w{z) = {4z + a ) {4z + h ) w{z) > (L4) 

and admits the series representation 

7-1 / t s (a)k{b)k z k , . 

2*1 M;c;*) = ^ , (1.5) 

fc=0 v 

where (a)^ = T(a + k)/T{a) is the Pochhammer symbol. 
The primary aim of this paper is to prove the following: 

• Theorem 1: 

The all-order e-expansions of the Gauss hypergeometric functions 

2-F1 (ii + as, J 2 + be; I 3 + ce; z) , (1.6a) 

2F1 (Ji + oe, I 2 + be; I 3 + \ + ce; z) , (1.6b) 

2-F1 (Ji + I + ae, I 2 + 6s; I3 + ce; z) , (1.6c) 

2 F 1 (Ji + \ + ae, I 2 + be; I 3 + \ + ce; z) , (1.6d) 

2F1 (Ji + 5 + ae, I 2 + 3 + oe; ^3 + | + ce; z) , (1.6e) 

where {Ik} are integer numbers, a,b,c are an arbitrary numbers, and e is an ar- 
bitrary small parameter, are expressible in terms of Remiddi- Vermaseren harmonic 
poly logarithms with rational coefficients. 

2. All-order e-expansion 

2.1 Non-zero values of the e-dependent part 

It is well known that any Gauss hypergeometric function may be expressed as a linear com- 
bination of two other hypergeometric functions with parameters differing from the original 
ones by an integer |26|, 27, p^, 25 1. Such a representation will be called a reduction, and 



the explicit algorithm will be called a reduction algorithm. Using the algorithm described 



in Ref. [25], the result of the reduction can be written as 



P(a,b,c,z) 2 Fi(a + h,b + I 2 ;c + h;z) = |Qi(a, b, c, z)-^- + Q 2 {a, b, c, z) j 2 Fi(a, b; c; z) , 

(2.1) 

where a,b,c, are any fixed numbers, P,Q\,Q 2 are polynomial in parameters a,b,c and 
argument z, and Ii,I 2 ,I 3 any integer numbers. 

All the hypergeometric functions ( |1.6| ) listed in Theorem 1 can be reduced to functions 
with I\,I 2 ,I 3 equal to zero for half-integer values of parameters, and to unity for integer 
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ones. In this way, all the hypergeometric functions of Theorem 1 are expressible in terms 
of the following five basic functions and their first derivatives: 



2 -Fi(ai£, a 2 e;l + ce;z), 2 Fi(aie, a 2 e; \ + fe; z), (2.2a) 
2 Fi(\+be,ae;l+ce;z), 2 Fi{\+be, ae; h+fe; z), 2 F l (\ + b l e,\+b 2 e;\ + f e; z). (2.2b) 



It was shown in Ref. fl2q| that only the two hypergeometric functions (and their first deriva- 
tive) of type ( |2.2a| ) are algebraically independent. The other three, (|2.2b| ), are algebraically 
expressible in terms of 2 Fi{a\e, a 2 e; \ + fe; z). Consequently, in order to prove Theorem 
1, it sufficient to show that the analytical coefficients of the e-expansion of the first two 
hypergeometric functions ( |2.2a ) are expressible in terms of Remiddi-Vermaseren polyloga- 



rithms. 

2.1.1 Integer values of of e-independent parameters 

Let us start from expansion of Gauss hypergeometric functions with integer values of 



parameters, and consider the function 2 Fi(a\e, a 2 e; 1 + ce;z). In ref. []30|, the all-order 



e-expansions for this functions and its first derivative were constructed in terms of multiple 



polylogarithms of one variable [16, 17]. These multiple polylogarithms may be expressed 
as iterated integrals 3 and have the expansion 



Li klM> ..., kn (z)= Yl fcl | fcT- ( 2 - 6 ) 

mi>m2>---m n >0 12 11 

Similar results were also derived (without explicit form of coefficients) in Ref. via 



nested sums approach. We will follow the idea of Ref. |30 |. 

The Gauss hypergeometric function 2 F\ (a±e, a 2 e; 1 + ce; z) is the solution of the differ- 
ential equation 

Tz ( z Tz + C£ ) w{z) = i z Tz + ai£ ) ( z Tz + a2£ ) w{z) ' (2 ' 7) 



3 Recall that multiple polylogarithms can be expressed as iterated integrals of the form 

f z dt dt dt dt dt dt dt dt 
Lifc, ... k„ (z) — / — o — o ■ ■ ■ o — o- o • ■ ■ o — o — o ■ ■ ■ o — o- , 2.3) 

u • n v ; J t t t 1 - t t t t 1 - 1 ' y ' 

k x -\ times k n -i times 

where, by definition 

f Z — o — o o dt o dt - f" dtl f 1 dt2 [ tk ~ 2 §^2-1 P" 1 ' 1 dtfc i (24) 

JO t t t 1 — t J tl J t2 J Jo 1 — 

k x -i times 



and u>i = -r 2 -, so that 

1 — 2 ' 



The integral (2. J ) is an iterated Chen integral (see also |32j]) w.r.t. the two differential forms u>q = dz/z 

i,— ,k n {z)= I 0Jq 1_1 wi • • -tJo" -1 ^! . (2.5) 
Jo 



- 4 - 



with boundary conditions w(0) = 1 and z-^w(z) 



z _ Q = 0. Eq. ( |2,7D is valid in each order 
of e, so that in terms of coefficients functions Wk{z) defined as 

oo 

w(z) =Y,Mz)e\ (2-8) 

k=0 

it can be written 

^ ~ Z ^lz (f^i) Wk ^ = ( ai + ° 2 ~ i) ^-^^ + a l a 2^fc-2(^) (2.9) 



for k > with 



u>£,.(z) = , k < 0. 
The boundary conditions for the coefficient functions are 



w k (0) = 



z4- z w k {z) 



z=0 



k>l, 
, fe > . 



Let us introduce a new function p(z) defined by 



Pi?) = Z ~j- W ( Z ) = Y2Pk( z )£ k > 



k=0 



where the coefficient functions satisfy 



d 



Pk{z) = z—w k {z) . 



The boundary conditions for these new functions follow from Eq. (2.11): 

Pfc(0) = , k>0. 
Eq. (|2.9j ) can be rewritten as a system of two first-order differential equations: 

(1 - Z )j^Pi( Z ) = ( 0l+a2_ ^) Pi-l( z )+ a l a 2Wi-2(z) , 

Z-^-Wi(z) = pi(z) . 
dz 

The solution of this system can be presented in an iterated form: 

dt 



(2.10a) 
(2.10b) 

(2.11a) 
(2.11b) 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



Pi(z) = (ai + a 2 -c) 
dt 



f z dt f z u„ 

/ z -pi-i(t)+aia 2 / ——Wi- 2 (t)-cwi-i(z), i>l, 

Jo 1 " Jo 1 ~ t 



Wi(z) 



t 



Pi(t) , i>l 



(2.16) 



4 We may note that 



2*1 



1 + a%e, l+a2£ 

2 + CE 



1 + ce \ - 

: 1 = 



Pk+2(z) 



aia 2 
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Taking into account that wq{z) = 1 and po(z) = (the e-expansion of p(z) begins with the 
term linear in e), we obtain the first few coefficients, 



Pl (z) = w x {z) = 0, (2.17a) 



a\a 2 
w 2 {z) 

aia 2 2 

w 3 (z) 



-ln(l-z) = H(l;z) , (2.17b) 

Li 2 (z) = H (0, 1; z) , (2.17c) 

7c i ln 2 (l -z)- cLi 2 (z) = lc H(l, 1; z) -cff (0, 1; z) , (2.17d) 

7c S 1 , 2 (z)-cLi 3 («)=7 c ir(0 > l,l;z)-cff(0,0,l;z) , (2.17e) 



aia 2 

where we have defined 7c = ai + a 2 — c, and Li n (z) and S a ^(z) are the classical and Nielsen 
poly logarithms [33, |3], respectively: 

-3o,6(«) = ^ Q _ jTj b , y d£ , S a ,i{z) = Ll a+ i (z) . 



The functions H(A;z) are the Remiddi-Vermaseren harmonic poly logarithms [20], and A 
is a multiple index including only entries and 1. 

From the representation ( 2,16| ) and result for the first few coefficients ( p. 17 ) we may 
derive the following observations: 

• Corollary 1: The all- order e-expansion of the function 2 i ? i(aie, a 2 e; 1 + ce; z) may 
be written in terms of harmonic polylogarithms Ht{z) only, where the multiple index 
A includes only the values and 1. 

• Corollary 2: The analytical coefficient of e k in the expansion of 2 i 7 i(ai£, a 2 e; 1 + 
ce; z) includes only functions of weight k with numerical coefficients. 

• Corollary 3: The non-constant terms of the e-expansion of 2 -Fi(ai£, a 2 e; 1 + ce; z) 
are proportional to the product aia 2 in any order of e. 

The first and the last statement follows from the representation 5 ( |2.16| ), the explicit value 
of coefficients functions Wk(z), k = 0,1,2 (Eqs. (|2.17|) ), and definition of harmonic poly- 
logarithms [2C[. The second statement follows from the form of the solution of Eq. ( 2,16| ). 

The relation between harmonic polylogarithms H^(z), with multiple index A including 
only and 1, and multiple polylogarithms of one variable (Eq. 2.6)) is well known [21] and 
follows from the proper definition (see Sec. 2 in Ref. p0| ): 

Li kl)k2 ,.., kn (z) = g( 0,0,-- - ,0, 1, 0,0, ••■ ,0, l,--- p,0,-j- ,0, 1] z) • (2.18) 
fci-i times fc 2 -i times fc„-i times 



3 The Corollary 3 follows also from general properties of hypergeometric functions. 
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By continued iterations of Eq. ( 2.16| ) and Eq. ( 2,17 ), we have reproduced all coefficients 



of the e-expansion of the Gauss hyper geometric function presented in Eq. (4.7) of [25]. For 
the coefficients functions Pi(z), u^z), and ps(z), we find a more compact form. We also 
obtain the higher-order terms u>^(z) and pe(z) of the e-expansion. The results are 6 



ai<i2 
a±a 2 



ln 3 (l-z) + (c7 c -aia 2 ) ln(l-z)Li 2 (z)+c 2 Li 3 (z) + (c^ c -2aia 2 ) 81,2(2;) , 

o 

c 2 Li 4 (z) - - (c 7c - a x a 2 ) [Li 2 (z)} 2 + 7 c 2 Si, 3 (z) + (c 7c - 2aia 2 ) S 2>2 (z) , (2.19) 



^ = ^ 3 M(l-z)-c 3 LU(z)-n% >3 (z) 
a\a 2 24 

-c(c7 c -aia 2 ) ln(l - z)Li 3 (z) -c (c^ c -2aia 2 ) S 2 , 2 (z) 



-7c (c 7 c-aia2) ln(l - z) 



ln(l-z)Li 2 (z) + Si, 2 (z) 



(2.20) 



w 5 (z) 
aia 2 



a\a 2 



7 3 S M (z) - c 3 Li 5 (z) - c 7 2 S 2 , 3 (z) - c(c7 c - 2aia 2 ) S 3j2 (z) 



+ (c7c - 01^2) 



7 c Li 2 (z) Si 2 {z) - -jcF^z) - cF 2 (z) 



(2.21) 



" 7 " 120 ln5(1 " Z) + c4Li5 W + C ^ S2 ' 3 ^ 



O 



ln 3 (l - z)Li 2 (z) + 31n 2 (l - z)S h2 (z) + 61n(l - z)Si, 3 (*) 



+ ( c 7c - oia 2 ) ln(l - z) 



7c ln(l-z)Li 3 (z) + 2cLi 4 (z) 



-(c- ai)(c- a 2 ) (c7 c - 2aia 2 )hi(l - z)S 2)2 (z) 



+aia 2 (c7 c - aia 2 ) 



~ ln(l - z) [Li 2 (z)] 2 - 2Li 2 (z) Si, 2 (z) + 2iq (z) 



+ (c7 c - 2aia 2 ) 



7 2 S M (z) + c 2 S 3j2 (z) 



where we have introduced two new functions: 



f' z dr 

F x (z) = / — ln 2 (l-x)Li 2 (x) . 
Jo % 
l' z dr 

F 2 (z) = / — ln(l - x)Li 3 (x) . 
Jo x 



(2.22) 

(2.23) 
(2.24) 



There is an algebraic relation 7 between these two functions: 



F 2 (l - z) = Fi{z) - 21nzSi, 3 (z) + 2S 2 , 3 (z) - Li 2 (z) Si, 2 (z) - lnzln(l - z)S 1<2 (z) 

11 1 

--ln 3 (l - z)ln 2 z - -lnzln 2 (l - z)Li 2 (z) + -C2hi 2 (l - z) \uz 

-C2Si, 2 (z)-C3Li 2 (l-z)-C5 , (2.25) 



6 The FORM[^ representation of these expressions can be extracted from Ref. pj[ 
7 We are indebted to A. Davydychev for this relation. 
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where 

Fi(l) = 2C3C2 - Cs ~ 2.9176809 

In this way, at the order of weight 5, one new function 8 , Fx, which is not expressible 
in terms of Nielsen poly logarithms, is generated by the Laurent expansion of a Gauss 
hypergeometric function with integer values of parameters. In general, the explicit form of 
this function is not uniquely determined, and the result may be presented in another form 
by using a different subset of harmonic polylogarithms. 

2.1.2 Half-integer values of of e-independent parameters 

Let us apply a similar analysis for the second basis hypergeometric function 



iFx 



axE, a 2 £ 



1 



+ fe 



(2.26) 



In this case, the differential equation has the form 



d 



,h vTz-\ + f£ ) w{z) = { z Tz + ai£ 



z— + a 2 e ) w(z) , 



with the same boundary conditions w(0) = 1 and z4^w(z) 



z=0 



(2.27) 

0. Using the e-expanded 



form of the solution, and noting that Eq. fl2.8|) , and in fact, Eq. ( 2.27 ) is valid at each order 
of the e-expansion, we may rewrite Eq. fl2.27|) as 



(1-z) 



d 



1 



dz 2z 



,ax+a 2 ) 



f 



z— ) Wi^i(z)+aia 2 Wi-2(z) .(2.28) 



Let us introduce the new variable y such that 9 , 
1 - 



y 



2-1 



(1-1/)' 



1 + ^ 4y 



l-z 



Ay 



dz 



i-y d_ 

l + y dy 



and define a set of a new functions Pi(y) so that 



10 



z±vh(z) 
and, as in the previous case, 



1-y d 



l + y dy 



y~r m(y) 



i-y 
1 + 2/ 



pi(y) 



(2.29) 



(2.30) 



p{y) = z -j- w ( z ) = Y2pk(y)t 



k=0 



Compare with results of [ [35| . 
9 The form of this variable follows from the analysis performed in Refs. 
10 We may note that 



jJFi 



1 + aie, l+a,2£ 
§+/= 



1 + 2/e 1 - y 
2z l + y 



E 



Pk+2(y) 

a\CL2 



(2.31) 
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In terms of the new variable y, Eq. ( [2.28 ) can be written as system of two first order 
differential equations: 



y-rPAy) = («i+«2 — — 
dy l + y 



Pl - 1 (y)+2f 



i - y l + y 



pi-i(y)+aia 2 Wi- 2 (y) 



y-i-w k (y) 

dy 



-Pk{y) 



(2.32) 



The solution of these differential equations for functions Wi(y) and pi{y) has the form 



Pi(y) 



Wi(y) 



dt 



2/^-2(01+02-/) ' 



1-t 



l + t 



Pi-i(t) - (oi+o 2 ) [wi-x(y)-Wi-\(l)} 



r ^ 

+0i0 2 J —Wi-2{t) 



i > 1 



y dt 



Pi(t) , i>l. 



(2.33) 



The point z = transforms to the point y = 1 under the transformation ( |2.29| ), so that 
the boundary conditions are 



w k (l) = 0, k>l, 
Pfe(l) = 0, k>0. 



(2.34) 



The first several coefficients of the e-expansion can be calculated quite easily by using 
w {y) = 1 and p {y) = 0: 



Pi{y) = wi{y) = 0, 
P2(y) 



a\a 2 

w 2 (y) 

a±a 2 



ln(y) = H(0;y) 



1 



-ln 2 (y) = -H(0,0;y) . 



(2.35a) 
(2.35b) 

(2.35c) 



Continuing these iterations, we may reproduce the coefficients of the e-expansion of the 
Gauss hypergeometric function ( p. 26 ) presented in Eq. (4.2) of Ref. p5|]. Since the length of 
the expressions obtained for the coefficient functions P3(y), W3(y), P4(y), ^(y), Ps(y), ^5(7/) 
is similar to those published in Eq. (4.1) of Ref. [25], we don't reproduce them here. 11 
The higher order terms of e-expansion are relatively lengthy and therefore will also not 
be presented here. Unfortunately, as in the previous case, we are unable to calculate the 
^-coefficient of e-expansion without knowledge of previous ones. 
From representation ( |2.33| ) we deduce the following result: 

• Corollary 4: The all-order e-expansion of function \2. 2(\ ) can be written in terms 
of harmonic polylogarithms H^{y) of variable y defined in l{2.29[) and multiple index 
A with entries taking values 0, 1 and —1. 



11 M.Y.K. thanks to M. Rogal for pointing out a mistake in Eq. (4.1) of Ref. 
coefficient should be "—2(3/ — a± — (12)" instead of "—2(3/ — 2ai — 202)"- 



In the e term, the 



-9- 



This statement follows from the representation ( 2.33j ), the values of coefficients functions 
Wk(z), k = 0, 1, 2 (see Eqs. ( 2.11D , ( |2.17| )), properties of harmonic poly logarithms, and the 
relation between powers of logarithms and harmonic polylogarithms. Also, Corollary 2 
and Corollary 3 are valid for the hypergeometric function ( 2.26j ). 

We would like to mention that, in contrast to the Eq. ( |2.16| ), Eq. ( [2.33 ) contains a 
new type of function, coming from the integral J f(t)dt/(l + t). Another difference is 
that the first nontrivial coefficient function, P2(y), is equal to ln(y), instead of ln(l — 
z), as it was in the previous case. It was shown in Ref. |20| that terms containing the 
logarithmic singularities can be explicitly factorised (see Eqs. (21)-(22) in Ref. f2(|), so 
that the coefficient functions, Wk(y) and Pk(y) from Eq. ( 2.33j ), have the form 



Li(j) (y)-Li (S) (1) 



w k {y) = ^c{s,a,k)ln k j (y) 
j=o 
fc-1 

Pk{y) = ^2c{s,a,k)ln k ~ J (y) Li/sx (y) - Li/*\ (1) 

3=0 



(i) 



(2.36) 



where c(s, a, k) and c(s, a, k) are numerical coefficients, s and a are multi-index, s = 
(si, • • • s n ) and a = (crj, • • • , a n ), belongs to the set of the square roots of unity, = ±1, 
and Li^ (y) is a coloured multiple polylogarithm of one variable [15, 16, 17], defined as 



Li 



E 



mi 



mi>m2>---m n >0 



nil m 



rrir. 



(2.37) 



It has an iterated integral representation w.r.t. three differential forms, 

0. 



dy 
= — , 

y 



a 



L0 n 



<?dy 
I - ay' 



a 



±1, 



so that, 
Li 



CT 1 ,CT 2 ,"- ,°k 

s l. s 2>'-' . 



(v) 



.31-1. 



S2-1 



L0 n 



Sfc-1 



cn<T2—<7k ' 



1 



(2.38) 



(2.39) 



3J. 



The values of coloured polylogarithms of unit argument were studied in Refs. 

2.2 Zero-values of the e-dependent part of upper parameters 

In the case when one of the upper parameter of the Gauss hypergeometric function is a 



positive integer, the result of the reduction has the simpler form (compare with Eq. Q2.1| )): 

P(b, c, z) 2 F 1 (I h b + h;c + J 3 ; z) = Qi(6, c, z) 2 iq(l, b; c; z) + Q 2 (b, c, z) , (2.40) 

where 6, c, are any fixed numbers, P, Q\, Q 2 are polynomial in parameters b, c and argument 
z, and I\,l2,h are any integers. 12 In this case, it is enough to consider the following two 



2 The proper algebraic relations for the reduction are given in Ref. 
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basis functions: 2-^i(l> 1 + ae; 2 + ce; z) and 2F1O-, 1 + ae; ^ + /e; 2). The e-expansion of 
this function can be derived from the proper solution given by Eq. (|2?nj ) or Eq. fl^ ), 
using the relations 



2F1 

and 
2F1 



1, l+a 2 e 
2+fe 



1 + ce d I a\e, d2£ 

irm 93-2-ri 1 , 

ai^o aia2£ az \ 1+ce 



1 + ce 



E 

fc=0 



1, l+a 2 e 

3 
2 



3 r + /e 



hm 5-3-2^1 

ai^O 2ai02£ dz 



aie, ci2e 
2 + / £ 



1 + 2/e 
2z 



E 

fc=0 



Pfc+2^) 




ai«2 


01=0 

( 


pk+2{y) 




a±a2 


ai=0 



(2.41) 



(2.42) 



where we have used the differential relation 



dz 



2F1 



a, b 



ab 



2F1 



1 + a, 1 + 6 
1 + c 



and the brackets mean that in the proper solution, we can put a\ = 0. The functions pk 
are given by Eq. ( |2.17| ) and Eq. ( |2.35 ), correspondingly. Due to Corollary 3, the limit 
a\ — > must exist. 

The case when both upper parameters are integers may be handled in a similar manner. 
Theorem 1 is thus proved. ■ 



3. Some particular cases 

3.1 The generalized log-sine functions and their generalization 



For the case < z < 1 the variable y defined in (2.29) belongs to a complex unit circle, y = 
exp(i0). In this case, the harmonic polylogarithms can be split into real and imaginary parts 
(see the discussion in Appendix A of Ref. ]2^| ), as in the case of classical polylogarithms. 



|3^ | Let us introduce the trigonometric parametrization z = sin 2 |. In this case, the solution 
of the proper differential equations ( |2.16| ) and ( |2.33| ) can be written in the form 

f e sin^ f e sin^ 

Pi(Q) = (ai+a 2 -c) dcf> |pi-i (<?!>) +aia 2 / d<f> 1 Wi-2 (<£)- cwi-i (9) , i>l, 

Jo cos § Jo COS § 

COS 

Wi (9) = I d<P—j-pi{4>) , i > 1 , (3.1) 



sin I 



and 



sin 2 f p cos 2 

Pi(0) = (ai+a 2 -/) / # |pi_i(^)-/ / # f-Pi-i (</>)+ a i fl 2 / d<M_ 2 (</>) , 

Jo cos 7f jo sin ^ jo 

^(0) = i>i, (3.2) 
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respectively. In the first case, the solutions of the system of equations (3.1) are harmonic 



polylogarithms with argument equal to sin 2 |. In the second case, the result contains the 
generalized log-sine functions |33|, 38, 39, 40 1 and some of their generalizations studied in 
Ref. (see also Ref. |4^]). For illustration, we will present a first several terms of the 
e-expansion 13 (see the proper relations, Table I of Appendix C in Ref. |23[|): 



2^1 



1 + aie, 1 + a 2 e 



sin 2 



;i + 2/e) 



sm ( 



x [9 + 2e{ 7/ (Ls 2 (ir-6)-9L e ) - f (Ls 2 (9)+9l b 



+e 2 { 2/ 72/ Ls 3 (0)+2 7/72/ Ls 3 (vr - 0)-/ 7/ Ls 3 (29) 



44/ 7/ [Ls 2 {9) L e -Ls 2 (tt-0) l e +9L e l e ] + 4/ 2 Ls 2 (0) Z e - 4 7 ;Ls 2 (tt - 0) L e 
1 

— ( 

6 



4-2/X + 2 7 2 0L 2 + ^aiM 3 4- 7 /72/<2 



+e 3 |^72/ [(ai + a 2 ) 7/ Ls 4 (vr - 9) + / 2 Ls 4 (0) - 3/ 7/ Lsc 2 , 3 (0)] - ^/ 2 7/ Ls 4 (29) 

+4 ai a 2 [2/Cl 4 (0) - 2 7/ Cl 4 (it -9) - /Cl 3 (9) 9 - 7/ Cl 3 (vr - 9) 9} 
+2 [/Z fl + v L e ] [/ 7/ Ls 3 (20) - 2 7/72/ Ls 3 (vr - 9) - 2/ 72/ Ls 3 (9)} 
+2 [/Z fl + lf L e ] 2 [2 72J Ls 2 (vr - 0) - /Ls 2 (29)} + aia 272 /0 2 Ls 2 (vr - 9) 

~a 1 a 2 fLs 2 (29) 9 2 -^a 1 a 2 9 3 [flg+ V Lg]--9 [/Z e + 7 /L ] 3 - 2 7/72/ vrC 2 [fl e +V L e] 



+a 1 a 2 (3a 1 + 3a 2 - 7f)9( 3 - 2(ai + a 2 ) 7/72/ vrC 3 J + C(e 4 ) 



(3.3) 



and 



/ aie,a 2 e 

2*1 1 



,2 



1 + aia 2 e 2 j ^9 2 



+e J 2/Ls 2 (0) #-2 7/ Ls 2 (vr - 0) #+4 7/ Cl 3 (vr 



)4-4/Cl 3 (0) + (3a 1 +3a 2 -7/)C 3 



+e 2 l 2 7/72/ Ls 3 (vr - 0) 0-/ 7/ Ls 3 (29) #+2/ 72/ Ls 3 (0) 0+ — a x a 2 9 A 



-2 [/Ls 2 W-7/LS2 (vr - 0)] 2 + 7 / 72 /^vrC 2 J + 0(e 3 



(3.4) 



where 



L e = In ( 2 cos - 1 , l e = In I 2 sin - ) , 



13 The FORM representation of these expressions can be extracted from 
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the generalized log-sine function is defined as 



Lsf } (<?) 



d(f) <f ln J ' 



-fc-i 



o 



2 sin — 
2 



Lsj (6) 



Lsf (6) 



(3.5) 



and we use the notation Lsc2 > 3(6 1 ) for the special combination (see Eq. (2.18) in Ref. j23|) 
Lsc 2 , 3 (6) = iLs 4 (26) - ±Ls 4 (6) + 2Ti 4 (tan §) - 2 In (tan f ) Ti 3 (tan §) 



In 2 (tan §) Ti 2 (tan §) - ±6* In 3 (tan f ) , 



(3.6) 



where the functions Tijv (z) are defined as 33] 



Ti N (z) = Im [Lijv (iz)] = f 



LiAr (k) - Li/v (-k) 



Ti/vr (z) 



— TljV- 



(s) 



(3.7) 

These functions receive special interest in physics through their role in the so-called 
"single-scale" diagrams, which depend only on one massive scale parameter. The massless 
propagator-type diagrams, bubble-type diagrams and propagator-type diagrams on mass 
shell all belong to this class. In particular, the single-scale diagrams with two massive 
particle cuts correspond to hypergeometric functions with value of argument equal to z = 
1/4. |E3, 22, U In this case, the value of the conformal variable y is equal to the primitive 
"sixth root of unity", y = exp (i-|). In contrast to the case in multiple polylogarithms 
(|2.6j ) of the primitive sixth root of unity studied in Ref. [44] and the more complicated 



case in coloured polylogarithms of the sixth root of unity studied by Broadhurst in Ref. 
1 16], the physically interesting case corresponds to coloured polylogarithms of square root 
( 2.39[) (harmonic polylogarithms) with argument equal to primitive sixth root of unity. In 



this case, some new transcendental constants, in addition to studied in Ref. [44| will be 
generated. The set of independent constants up to weight 5 was constructed in Refs. 
|3l 



,0- 



3.2 Special cases: all-order e-expansion in terms of Nielsen polylogarithms 

One advantage of a trigonometric representation used in the previous section is the theo- 
rem proved in Ref. || (see also Ref. |4(|), that any generalized log-sine function (|3.5|) is 



expressible in term of Nielsen polylogarithms [34] only. Using this theorem, it was shown 
in Ref. 



23] that for the Gauss hypergeometric function 



2F1 



1, 1+ae 



sin 



2 6 



(3.8) 



the Laurent expansion is expressible in terms of only Nielsen polylogarithms in the three 
cases (i) 6 = 0, (ii) b = a, (hi) a = 2b. Using the reduction algorithm [25], we can claim that 
the Laurent expansions of the following functions are also expressible in terms of Nielsen 
polylogarithms only: 



2F1 



h,h+e 



sin ^ 



h,h+£ 



sin 5 



2F1 



h,h+£ 



sin 



,(3. 
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where ii, h and I3 are integers. 

It is interesting to analyze this solution from the point of view of Eq. ( |3.2[) . Due to 
fact that a\ = 0, the last term in Eq. ( |3.2|) is identically equal to zero. In case (i), only 
the first term survives, with integration kernel having the form dln(cos^). In case (ii), 
only the second term survives, and the integration kernel has the form <iln(sin ^). In case 
(hi), the first and second terms can be reduced to the second case of a double argument. 
The statement about expressibility of inverse binomial sums in terms of log-sine function, 
proved in Ref. (2^] (see also Ref. [^|]) applies to all three of these cases. 

We can extend the class of Gauss functions whose e-expansions are expressible in terms 
of only Nielsen polylogarithms by using algebraic relations 14 between of the fractional-linear 
arguments (see Sec. 3 in Ref. p5| ). The cases which may be expressed in this manner are 
summarized in Table I, where a, b, c are parameters of the Gauss hypergeometric functions 
2-Fi(a, b; c; z) and h, h and ^3 are integer: 



Table I 


a 


h 


h 


h 


h 


h 


h 


h+e 


h+e 




b 




\+h+e 


l+h-e 


\+h+e 


\+h 


\+h+e 


\+h 


\+h+e 


\+h+e 


c 


\+h+e 




\+h+e 


h+e 


h+e 


h + 2e 


h + l+e 


h + l+e 


h + l + 2e 



The results of this section can be formulated as follows: 



• Proposition 1: All cases of Gauss hypergeometric functions with half-integer val- 
ues of parameters for which the all-order e-expansion is expressible in terms of only 
Nielsen polylogarithms are described in Eq. \3. !\ ) or the parameters shown in Table I. 



4. Conclusions 



The main result of this paper is the proof of Theorem 1, as stated also in the abstract. 
The proof includes two steps: (i) the algebraic reduction of Gauss hypergeometric functions 
of the type in Theorem 1 to basic functions and (ii) the iterative algorithms for calculating 
the analytical coefficients of the e-expansion of basic hypergeometric functions. 

In implementing step (i), the algebraic relations between basis functions with half- 
integer values of parameters reduce all of the cases to the one basic function of type ( 2.2a| ) 



and its first derivative (see details in Ref. [^] ) . In step (ii) , the algorithm is constructed for 
integer values of parameters in Eq. ( |2.16 ) and for basis Gauss hypergeometric functions with 



half- integer values of parameters in Eq. (|2.33[) . This allows us to calculate the coefficients 



directly, without reference to multiple sums. 

It is interesting to note that the Laurent expansions of the Gauss hypergeometric func- 
tions with integer values of parameters are expressible in terms of multiple polylogarithms 
of one variable (see Eq. fl2.6|) ) or the Remiddi-Vermaseren harmonic polylogarithms with 
multiple index including only values and 1 . The argument of the resulting functions coin- 
cides with the original variable of the hypergeometric function. For Gauss hypergeometric 



4 This can also be derived via the integral representation. 
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functions with half-integer values of parameters, the coefficients of the e-expansion produce 
the full set of harmonic polylogarithms, or coloured multiple polylogarithms of one variable 
(see Eq. ( 2.37Q ). These functions depend on a new variable, related to the original variable 
by conformal transformation (see Ref. [25]). 

For special values of the argument of the hyper geometric function, z < 1, the coloured 
multiple polylogarithms of one variable may be split into real and imaginary parts. This 
case has been discussed in section 3.1. It was shown that the physically interesting case, 
representing single-scale diagrams with with two massive particle cuts, corresponds to 
coloured polylogarithms ( |2.39|) with argument equal to a primitive "sixth root of unity" , 
y = exp (i^)- This gives an explanation of the proper "basis of transcendental constants" 
constructed in Refs. [38] and ||, and its difference from the proper basis of David Broad- 
hurst [fill. 



In the section 3.2, the subset of Gauss hyper geometric function is analyzed, showing 
that the all-order e-expansion is expressible in terms of Nielsen polylogarithms only. In par- 
ticular, we have formulated the proposition that the only Gauss hypergeometric functions 
with half-integer values of parameters for which the all-order e-expansion is expressible in 
terms of Nielsen polylogarithms only belong to one of the functions described in (3.9) or 
in Table I. 

In Appendix we discuss the construction of the all-order Laurent expansion of the 
Gauss hypergeometric function fl2.26j ) around z = 1. 



Acknowledgments 

We are grateful A. Davydychev for useful discussion. M.Yu.K. is thankful to participants of 
conference "Motives and Periods", University of British Columbia, Vancouver, June 5-12, 
2006 [j49| , for interesting discussions. Special thanks to Andreas Rosenschon for invitation 
and financial support and D. Kreimer and H. Gangl for enormously and useful discussions 
and suggestions. M.Yu.K. is very grateful to Laura Dolchini for moral support when paper 
was written. This research was supported in part by RFBR grant # 04-02-17192 , NATO 
Grant PST.CLG.980342 and DOE grant DE-FG02-05ER41399. 

Appendix 

A. The Laurent expansion of Gauss hypergeometric functions with half- 
integer values of parameters around z = 1 

The identities between harmonic polylogarithms ( |2.6[) under the action of the group of 
fractional- linear transformation of the argument, 

z — > 1 — z, -, z 2 , 

z 

was considered in Ref. |2(J (see also Ref. [^]). It was shown |2^] that the full set of 
Remiddi-Vermaseren functions is invariant with respect to transformations 

1 1 - z 

T+z~- 



- 15 - 



In this Appendix, we wish to show that the coefficient functions uJk(y) an d pk(y) entering in 
the Laurent expansion of the hyper geometric function ( p. 26 ) satisfy to some identities with 
respect to the argument transformation z — > 1 — z. The derivation is trivial if we recall that 
the all-order Laurent expansion of the hypergeometric function ( 2.26j ) can be written in 
terms of coloured polylogarithms of argument y. Under the transformation z — » 1 — z, the 
conformal variable simple changes its sign, y =— ► —y. The full set of Remiddi-Vermaseren 
functions is invariant (up to the addition of a constant imaginary part) with respect to 
changing of the sign of this variable. Consequently, the coefficient functions Wfc(y) and 
Pk(y) should be related to u>k{—y) and Pk(—y)- 

Let us present the explicit relations. Using the Rummer's relations between hyperge- 
ometric functions of arguments z and 1 — z, we obtain 



2^1 



/ a\E,a2E 


- 







z l/2+(f-ai-a2)e T (|+/g) £ Q + ( Ql + q 2 - f)e) d 

(l-2)-V2+/e r(l+aie)r(l+a 2 e) dz 

r(|+/ £ )r (|+(/- a x -a 2 )e) 



2F1 



— a\E, —(l2£ 

i + (/-ai-a 2 )e 



(A.l) 



T (!+(/ - oi)e) T (| + (/ - a 2 )e) Zi 1 ^ i + (ai + a 2 - f)e J ' 
The all-order e-expansion for the hypergeometric functions entering in r.h.s. of this relation 



is constructed in Sec. |2.1.2j . Let us apply the same technique for constructing the Laurent 
expansion of the hypergeometric function on the l.h.s. In accordance with standard proce- 
dure (see Ref. [pq] ), let us introduce a new variable, Z = 1 — z, so that the the differential 
equation around Z = has the form 



This equation is equivalent to Eq. ( |2.27| ) with the proper change of variable and one of the 
parameters, 

(z,f) < — * (Z,ai +a 2 - f) , 



so that we can use the results of Sec. 2.1.2| with the proper change of notations. In partic- 



ular, the solutions of the differential equations for the functions Pi(Z) and Wi(Z) have the 
form 

1 . 1 



Pi(Y) 



Y 



dt 



+a\a2 



-Y 



2(01+03-/) 



l-t 



-V 



1 + t 



Pi-i{t) - (01+02) [wi-i(Y)-Wi-i(l)] 



dt 

T 



Wi-2(t) , i > 1 , 



wi(y) = -] jPiV), *>i- 

where new variable V is defined as 



g 

z-i 



1 + 



z 

Z-l 



(A.3) 



(A.4) 
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and y is defined by Eq. ( 2.29 ). In this way, both parts of relation ( |A.1| ) are expressible in 
terms of coloured polylogarithms depending on the arguments y (r.h.s.) and — y (l.h.s.). 

We expect that relations following from Eq. ( |A.1| ) may be useful in obtaining some dual 
relations for coloured polylogarithms (see Ref. (if], |l7|, Q), and in the obtaining algebraic 
relations between coloured polylogarithms of the primitive "sixth root of unity" , as in case 
of multiple zeta-values [47]. At the present moment, we are not ready to discuss these 
relations. 
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